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Abstract. The quantized electromagnètic field inside and outside an absorbing 
high-Q cavity is studied, with special emphasis on the absorption losses in the coupling 
mirror and their innuence on the outgoing field. Generalized operator input-output 
relations are derived, which are used to calculate the Wigner function of the outgoing 
field. To illustrate the theory, the preparation of the outgoing field in a Schròdinger 
cat-like state is discussed. 
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1. Introduction 

Cavity QED provides a rich tool for experimental tests and applications of quantum 
mechanics, with special emphasis on quantum decoherence and communication [1]. The 
use of high-Q cavities makes it possible to control the evolution of coupled atom-field 
systems in order to synthesize, at least on principle, arbitrary nonclassical states of light 
[2, 3, 4]. For example, schemes have been developed to generate various nonclassical 
states of light by mapping atòmic Zeeman states on states of cavity fields [8, 9]. 
Nonclassical states of light on their part play a fundamental role in quantum information 
processing [5, 6, 7]. 

Unfortunately Nature imposes limits on us, because of quantum decoherence, which 
is unavoidably connected with the losses that are always present in practice. In cavities 
both wanted and unwanted losses play a role. Whereas the wanted losses, which result 
from the input-output coupling, typically affect the quantum state of the cavity field, 
unwanted losses such as scattering and absorption losses affect both the quantum state 
of the cavity field and the quantum state of the outgoing field. Keeping track of leaking 
emission renders possible for real-time measurement of the quantum state of the field 
and, therefore, offers the possibility of feedback control of open cavities [10], provided 
that the input-output coupling can be regarded as being the main source of decoherence. 
For high-Q cavities this must not necessarily be the case. Moreover, when in quantum 
networks the radiation emitted by cavities plays the role of quantum communication 
channels connecting, e.g., àtoms trapped in distantly separated cavities [11, 12], even 
small unwanted losses can lead to a noticeable degradation of the nonclassical features 
of the radiation [13]. 

In the rest of this article we study the problem of the absorption losses of a high-Q 
cavity and their effect on the quantum state of the outgoing field in more detail, with 
special emphasis on the absorption losses in the coupling mirror. Basing the calculations 
on quantum electrodynamics in dispersing and absorbing media, we derive a formula 
that relates the Wigner function of the outgoing field to the Wigner functions of the 
cavity field, the incoming field, and the dissipative channels. As an application of the 
theory, we study a scheme for generating Schròdinger cat-like states by combining two 
incoming modes prepared in coherent states with a cavity mode prepared in a squeezed 
number state. 

The paper is organized as follows. In section 2 the cavity model is described and 
the bàsic equations are given. Expressions for the cavity field and the outgoing field are 
given in sections 3 and 4, respectively. The Wigner function of the quantum state of 
the outgoing field is derived in section 5, and section 6 presents an application. Finally, 
the main results are summarized in section 7. 
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Figure 1. Scheme of the cavity. The fractionally transparent mirror (region (2)) 
is modelled by a dielectric plate. The active sources are in region 1, which can also 
contain some medium. 



2. Cavity model 

For simplicity let us consider a modified version of Ley and Loudon's [14] cavity model, 
defined by a planar multi- layer system according to figure 1, where the layer labelled 
j = is assumed to be the totally reflecting mirror and the semitransparent mirror 
(layer 2) is modelled by a dielectric plate. The linearly polarized electromagnètic waves 
propagate in the z direction, for which we use shifted coordinate systems such that 
< z < l for j = 1, 0<z<d for j = 2, and < z < oc for j = 3. Applying the one- 
dimensional version of the quantization scheme in reference [15, 16], allowing for active 
atòmic sources inside the cavity [Ath atom with elèctric dipole moment aU(í) being at 
position za\, and following reference [17], we may represent the elèctric field in the jth 
layer of permittivity Sj{uj) (j = 1, 2, 3) in the form of 



/■oo 

È ü \z,t) = / dujÈ {j) (z,u,t) + H.c, 
Jo 

M j) (z, cu, t) = ÈÜl(z, u, t) + ÈÜ\z } 00, t), 



where 



•0'')/ / 



È£ J ee (z,u,t) = iufio I dz' G ( "'\z,z',u) j^J e (z',u,t) 



(1) 
(2) 

(3) 



V / dz'G^' 

is the free- field part indicates integration over the j'th layer), 

2 c 

È[ j \z,u,t) = ^-^-Y / dí'eíí-íV^'"* ^') lmG^\z A ,z,Lü){4) 
ne A c 2 



(A: mirror area) is the source-field part, and 



He, 



TïA 



(5) 
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[£ ) e (z,^í) = e--( í - í ')/f2 c (^^>0; f U) (z,u,t), bosonic fields; e» = e» + ie»]. 
Further, 

z', cu) = ^[£ j> (<z, cu) cu)6(j - /) 

+ ^<(^c)S^">(^,cu)0(j"-j)] (6) 
is the (nonlocal part of the) Green function, where the functions 

£<i>(z, lo) = é^-^ + r j/3 e' l/3 ^ z - d ^ (7) 

and 

£ ü)< (z, cu) = e~^ z + r j/0 e^ z (8) 

represent waves of unit strength travelling, respectively, rightward and leftward in the 
jth layer and being reflected at the boundary [note that Q(j — f) means Q(z-z') for 
j=j'}- Further, E"' is defined by 

E jf = 1 ton* 3 3 *3/i' e J 3 ? (g) 
Dj = l- r j/0 r j/3 e 2i ^ (10) 



where 



and 

Pi = PM = v^M 7 = + 7 = /?■ + (n) 

> 0, di = Z, d 2 = d, d 3 = 0]. The quantities tj/j' = (Pj/ and r^/ are 
respectively the transmission and reflection coefficients between the layers f and j, 
which can be recursively determined. Note, for the cavity model under consideration 
río — — 1, since perfect reflection from the left-side mirror has been assumed. 



3. Cavity field 

To evaluate, according to equations (l)-(4) for j = 1, the elèctric field inside the cavity, 
we notice that the function Di(u) given by equation (10) defines the spectral response 
of the cavity. That is, the complex resonance freqüències fl k = u k — \i^k of the cavity 
field are determined by the zeroes of 

D 1 (Ü k ) = l + r 13 (n k )e 2 ^^ l = (12) 

and may be found by iteration. Somewhat lengthy but straightforward calculations 
show that - provided that T k <C Au k , with Au k = \{uj k+ i — cü k -i) being the width of 
the fcth interval - equation (1) [together with equations (2)-(4)] can be approximately 
rewritten as 

(z, t) = J2 E k (z)à k (t) + H.c, (13) 
k 
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E k (z) = iu k 



sin[/3i(w fc )z], 



and 



€oei(üJ k )lAu) k _ 

T k b kin (t') + J2 A k^kx(t') 



(14) 



,2ni(cj fe )/ 

Here the operators c^í) and b k i n (t) are defined by 

1 



^£ fc (z A )ík(0}.(15) 



fefcin(í) = 



7 27T J(A fc ) 

1 /■ 



dí^ c A (w,í), 
du;ò in (u;,í), 



(16) 
(17) 



'2-K J {Ak ) 

where the integration runs in the interval (A k ) = [^(uj k -i + cu k ), \{oj k -\-üu k+ i)], and the 
bosonic operators 6i n (u;,t) and è\(u,t) (A = cav, ±) read 



b in (u,t) = - 



HocttA \n 3 



dz'e^'^Jz',u;,t), 



c±(uj, t) = —at± 



flQCTlA 1 



^ 2n2 J[2] 



àz' 



~ e -i/3 2 (z'-d) ± e if3 2 z> 



^Jz',u;,t), (19) 



C C£LV (UJ, t) = - a cav\/^ ~7^ y dz SÍn (Pl Z )íhL( Z > U > *)' 



'[1] 



with 



«cav = «cav(^) = 2«s/2|m| K sinh(2^'/) - < sin(2#Z)] " ; 
a± = a ± (w) = \n 2 \e^ d/2 [rí 2 sinh((3 2 d) ± n 2 sm((3' 2 d)f^ , 
and the coeficients T k , A k \ are defined according to 



T k = T(üü k ), T(u) = 

A k ± = A±(u; k ), A±(u) = 



t2i\/nï 



D' 2 a ± 



(r 23 e^ d ± 1) 



A kcw = A cav (uj k ), A cav (uj) = -Ai 



Oir 



(20) 

(21) 
(22) 

(23) 
(24) 
(25) 



It can be proved that the operators b kin (t) and c k \{t) satisfy, on a (course-grained) 
timescale Aí^> AccÇ 1 , Aü;^/ 1 , the commutation relations 

ihUt),bljt')] = 5 kk ,5(t-t'), (26) 

[èkx(t), cl, x ,(t')\ = 5 kk ,5 xx ,5(t - t'), (27) 
and the equal-time commutation relation 

[à k (t),àl(t)]=8 kk , (28) 
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hold. From equation (15) it follows that à k (t) obeys the Langevin equation 



i 

c 1 2 r_ í 

+ 



à k = -i(oj k - \iY k ) à k - - ^ E k (z A )d A 



2ni(u k )l 



T k b ki n(t)+J2 A kx è kx(t)\. (29) 

A 

Within the approximation scheme used, the total decay rate T k determined from 
equation (12) can be decomposed into a radiative part 7fc ra d and a non-radiative part 

7fcabs i 

Tfc = 7fcrad + 7fcabs, (30) 

where 

7fcrad= 2K(L fc )|/ |Tfc|2 ' (31) 

7, abs = $>a = ^J^i Ç l^l 2 . (32) 

Note, that Langevin equations of the type of equation (29) can be also obtained from 
quantum noise theories on the basis of an appropriately chosen Senitzky-Gardiner- 
Collett Hamiltonian [18, 19]. 

4. Input-output relations 

For simplicity we restrict ourselves to the case of a cavity in free space, i.e., n 3 — > 1. 
To determine the input-output relations, we decompose the field outside the cavity 
[equations(l)-(4) for j — 3 together with the Green function (6)] into incoming and 
outgoing parts and introduce the bosonic operators 

/ TïA - (3) 

b out (uj,t) = 2^—^E OVit (z,uj 1 t)\ z=Ç)+ . (33) 

Recalling the resonance structure of the spectral response of the cavity, we may write 

6 out (í) = — ]= i dujb out (uj,t) = y2b kout (t), (34) 
V2tt J ^ 

6fcout(*) = — ij== I dub kout (u,t), (35) 



'2TÏ J( Ak ) 

where, within the approximation scheme used, 



b k0 ut(u, t) — 



nA a; fc íi3exp[i/3i(íí fc )Z] 



x 

A 



fxochuj e ei(Q k )lA 
í àt'Q{t-t')e- iü ^d A {t')sm[mk)z A \ 

4 J 

+ 2^ T(0) H / d^e(í-Oe^ (t -°[T(ü;)Uw,0+Z)^H 6 >í ü; ' í ')' 

À 

+4 o) (ü;)c + (w ) í) + #(w)c_(w, í) + i?^(a;)6 in (a;, í), (36) 
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with 

^± M = — , 37) 

D 2 a± 

T W( w ) = ÍlLe lAÍ , (38) 
v n i 

R {o) (uj)=r 31 . (39) 

Performing in equation (35) the lü integration and recalling equations (15), (16), 
and (17), we obtain the generalized input-output relations 



'fcout 



(*) 



^-^j TÍ°)a fe (0+4 o) ^n(í)+4 O -íc fe+ (í)+4 O 2c fe -(0 (40) 

[T fe (o) =TV>(u k ), A^l = A^\u k ), R k ° ] = R io) (tü k )]. In close analogy to equation (26), the 
commutation relation 

lL·ut(t),bl out (t')} = S kk ,5(t-t') (41) 

can be shown to hold. Omitting the last two terms on the right-hand side in equation 
(40), which describe the effect on the outgoing field of the noise associated with the losses 
in the coupling mirror, we recover the input-output relations suggested by Standard 
quantum noise theories [19]. Note, that a Hamiltonian of the Senitzky-Gardiner-Collett 
type [18, 19] is not suited for consideration of the effect. 

5. Quantum- state extraction 

Let us suppose that at some initial time ío the kth cavity mode is prepared in a certain 
quantum state and evolves freely for t > to (the preparation time is assumed to be short 
compared to the decay time T^ 1 of the cavity mode). From equations (29), (36), and 
(40) it then follows that 

fytout(^, t) = F*(lü, t)à k (t ) + B k (cü, t), (42) 



where 



V2 

1 I C \ ~-,(o)* 



with 



^^vstej T > M ^ (43) 

È k (u, t)= du' \G* kin (u>, J, t)b kin (u', í )+V G* kX (uj, J, t)c kX {u', t )] ,(44) 



Lü — \i k 

v{u,u',t) = Í--f--J_[/(a;,í) - f(u',t)e^'-^}, (46) 

Z7T ln\l Lü — Lü 

G kin (u, u', t) = T^T^u, lü', t) + R^*e^-^S{u - lü'), (47) 

G kcav (u, lü', t) = T^*A* kcElv v k (u, lü', t), (48) 

G k ±(cü, lü', t) = T^*Al ± v k (Lü, lü', t) + A k fe^^- to) 5(Lü - lü'). (49) 
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For the following it will be convenient to introduce unitary, explicitly time-dependent 
transformations according to 



i 

c kX (u,to) =^2<f)^*{uj,t )c ( il{t ), 



(51) 
(52) 



where, for chosen t and ç (ç = out, in, À), the functions (f>ç\u,t) are complete sets of 
square integrable orthonormal functions. Needless to say, that the bosonic commutation 
relations are preserved. 

From equation (42) it is easy to see, that the cavity mode predominantly escapes 
into the outgoing mode with 



Vk{t) 



where 



Vk(t) 



du\F k (u,t)\' 



(A fc ) 



1/2 



Hence, for chosen k, equation (42) takes the form 



where 



feout 



Bf{t) 



(t) 



Vk(t)à k (t ) + È®{t) ifi = l, 
B^\t) otherwise, 

d ^ <f>olt(u,t)B k (uj,t), 



(A fe ) 



(53) 
(54) 

(55) 

(56) 

(57) 
(58) 



and it can be proved that the commutation relations 

[à k (t ),è™(t)]=Q (f >í„), 

[B£\t),èW\t)] = l-ril(t) 
hold. 

It is now straightforward to calculate the characteristic function C^ ut (/3,t) of the 
quantum state of the relevant outgoing mode; for details, see reference [20]. The result 
reads [a = (in, A)] 

C2L(M = exp[-I|/5| 2 [l - V *(t) - \xïl(t)\ 2 { 

<T,Í 

x c^mXlctlixtm, (59) 

ai 

where the mode couplings are defined as 

x£(*)= í [ du'G ka (u,u f ,t)<f>W\u,t)<f>?(u;',t ), (60) 

J(A k ) J(A k ) 
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Cfc(/3) is the characteristic function of the initially prepared quantum state of the cavity 
mode, and C^((3) are the characteristic functions of the modes of the incoming field 
(o" = in) and the dissipative channels (er = À). From equation (59) the corresponding 
relation between the Wigner functions is derived to be [^(í) = r] k (t), x^\t) = Xkl(t)] 



2 I 

W out (a,t) = -- 



x 



j à 2 a'W{a')\[ J d 2 a ai W^(a ai ) exp 



,(6i) 



where (7 rad = 7 fcrad , 7abs = 7fcabs) 

(o) 

fj 2 ^ — ^ad ^ _ e -(7rad+7ab S )(í + Aí-í(,)j ^ ^) 

7rad "I - 7abs 

with 

,>)_,>) _ C i r (°)|2 / fi o\ 
7rad — 'ferad ~ 2\m\l k 

In particular, when the incoming modes and the dissipative channels are in thermal 
states (with average number of thermal photons n ai ), equation (61) reduces to 

2 1 

W out (a,t) = — 



x J d 2 a' exp 



Tl-^(t) + 2D <ri n a< |xí? ) (í)l 2 
2|77(í) a' -a\ 2 



i-77 2 (0 + 2E^lx^(0l 2 

which reveals that the condition 

^ n— » i (65) 

must be satisfied to ensure that the quantum state of the outgoing field becomes 
sufnciently close to the quantum state that the cavity field was initially prepared in. 
It should be pointed out that the additional noise terms in the input-output relation 
(40), which are associated with the losses in the coupling mirror, modify the functions 
Gka(^, w', t) [equations (47)-(49)], which leads, compared to Standard noise theory, to an 
increased value of 2 n ci \xa\t)\ 2 in equation (64), with the result that the condition 
of almost perfect quantum state extraction becomes more restrictive. 

6. Generation of Schròdinger cat-like states 

As an application of the theory, let us consider the case where two incoming modes 
with Xü?(t) = Xin{t) = Xm(í) are initially prepared in coherent states \f3) and \—i/3), 
respectively, and all other input modes and the dissipative channels are in thermal states 
with average numbers of thermal quanta n. Then, from equation (61) we find 

2 1 

W out (a,t) = — 

7r 1 — r] z (t) + 2n\x{t)\ 



W(a'), (64) 
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Figure 2. Wigner function of the output state for n = 1, r = 0.8, n — 0.001, /3 = 2, 
(a) x = 0.1, Xin = 0.9, r] = 0.9, (b) X = 0.7, Xin = 0.7, 77 = 0.7. 

where the notation 

Y^^\t)\ 2 = n\ X it)\ 2 (67) 

ai 

has been introduced, with the two incoming modes being excluded in the primed sum. 
Let us suppose that the cavity mode is initially prepared in a squeezed number state \r, n) 
= S{r)\n) [S(r), squeeze operator]. Examples of the Wigner function of the quantum 
state of the outgoing mode are plotted in figures 2 and 3 for r = 0.8, n = 1 and 
r = 1.1, n = 10, respectively, with n = 0.001, (3 = 2 for both cases. As we can see 
from figures 2(a) and 3(a), the quantum state of the outgoing mode exhibits the typical 
features of a Schròdinger cat-like state, provided that the mode coupling is strong enough 
(Xin - > 1) and the noise associated with the losses is sufliciently small (77 — > 1, x — > 0). 
Figures 2(b) and 3(b) reveal, that with decreasing strength of the mode coupling and/or 
increasing noise the nonclassical features of the quantum state of the outgoing mode are 
completely lost. For the parameters chosen the nonclassical oscillatory fringes typical 
for a Schròdinger cat-like state can be observed as long as rj(t) > 0.7, which for t — > 00 
corresponds to the requirement that 7 a b s ^ 7rad- Note that the additional noise terms 
in the input-output relation (40), which are associated with the losses in the coupling 
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Figure 3. Wigner function of the output state for n = 10, r = 1.1, n = 0.001, (3 = 2, 
(a) x = 0.1, Xin = 0.9, rj = 0.9, (b) X = 0.7, Xin = 0.7, V = 0.7. 

mirror, effectively reduce the strength of the mode coupling Xm- However, the necessary 
mode analysis will be considered in detail in a forthcoming paper. 

7. Summary 

To summarize, we have studied the input-output problem for a high-Q cavity, with 
special emphasis on used input channels, taking into account the absorption losses in 
the coupling mirror. Within the framework of exact quantization of the electromagnètic 
field in dispersing and absorbing media we have performed the calculations for a one- 
dimensional cavity bounded by a perfectly reflecting mirror and a semi-transparent 
mirror. The theory generalizes the Standard quantum noise theory based on a 
Hamiltonian of Senitzky-Gardiner-Collett type [18, 19], which does not fully take 
into account the absorption losses in the coupling mirror. Using the generalized 
operator input-output relations, we have calculated input-output relations for the 
Wigner functions of the quantum states involved. Finally, we have applied the theory 
to the case where two incoming modes prepared in coherent states are combined with a 
cavity mode initially prepared in a squeezed number state to prepare the outgoing field 
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in a Schròdinger cat-like state. 
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